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Abstract: We systematically study the spectrum of open strings attached to half 
BPS giant gravitons in the N = 4 SYM AdS/CFT setup. We find that some null 
trajectories along the giant graviton are actually null geodesies of AdS^ x S 5 , so that 
we can study the problem in a plane wave limit setup. We also find the description 
of these states at weak 't Hooft coupling in the dual CFT. We show how the dual 
description is given by an open spin chain with variable number of sites. We analyze 
this system in detail and find numerical evidence for integrability. We also discover 
an interesting instability of long open strings in Ramond-Ramond backgrounds that 
is characterized by having a continuum spectrum of the string, which is separated 
from the ground state by a gap. This instability arises from accelerating the D- 
brane on which the strings end via the Ramond-Ramond field. From the integrable 
spin chain point of view, this instability prevents us from formulating the integrable 
structure in terms of a Bethe Ansatz construction. 
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1. Introduction 

The gauge theory /gravity duality is probably one of the deepest ideas in theoretical 
physics, in that it gives us in principle the possibility to understand quantum gravity 
exactly in some setups. The main tool to address this duality involves the ideas of 't 
Hooft of the large N expansion of field theories in terms of a dual string theory [p]] . 
A lot of recent progress has happened because there are some examples where the 
dual string theory is known. The simplest and most studied example of this duality 
is the original formulation of the AdS/CFT correspondence in terms of maximally 
super symmetric Yang Mills theory in four dimensions and type IIB superstring theory 
compactified on AdS 5 x 5 5 ||. 

The setup is such that in the limit where R (the radius of curvature of AdS§ x 5 5 ) 
is large, the field theory is strongly coupled in the sense of 't Hooft: i? 4 ~ gy M N 
is large. Indeed, it would be very nice to have a good understanding of how this 
duality works in detail, as it would also give us hints on how to calculate observables 
in strongly coupled gauge theories in four dimensions for more general setups. 

There are various avenues to explore this correspondence, based on different 
observables that one would want to work with. Historically, the first set of observables 
to be addressed were the dual states to single graviton perturbations of AdS 5 x 5 5 
[0, M and how to go about testing their correlation functions. These states are all 
members of BPS multiplets, so their energies are protected by supersymmetry. 

A few years later, it was discovered that there were interesting geometric limits 
where the string theory could be quantized exactly |§, |§, 01- Translating the cor- 
responding quantum numbers of states to the field theory language provided a new 
large N limit where one also scales the energy and R charge of the observables as one 
makes the coupling constant large ||. The subsector of states that one focuses on is 
closely tied to the supersymmetry of the original system, so the states in questions 
are nearly supersymmetric. In these cases, perturbation theory at strong coupling 
was parametrically suppressed by the large quantum numbers of the states in ques- 
tion, so it was possible to use perturbation theory to examine the strong coupling 
regime of the field theory in a relatively safe environment. 

This result produced an interest to make a systematic study of perturbation 
theory near the free field limit. In particular, it became interesting to compute the 
full spectrum of anomalous dimensions of all local operators on the field theory. 

In the case where one focuses on operators whose free field dimension grows 
at most as yN, the states are roughly described by a Fock space of closed string 
states, and the dimension of the operators in the large N limit is dominated by 
the planar diagram expansion. The expansion to one loop order of this problem 
revealed a very surprising structure. Minahan and Zarembo discovered that the 
one loop spectrum of anomalous dimensions in a subsector of the theory gave rise 
to an integrable 50(6) spin chain: a generalization of the Heisenberg XXX spin 
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chain. This result was later generalized to include the full set of local operators of 
the theory and it was shown that the full planar one loop spectrum of anomalous 
dimensions gave rise to an integrable spin chain model that could be solved by Bethe 
Ansatz techniques [ITJ. In a parallel development, Bena et. al. [1 1] discovered that 



the string sigma model on AdS$ x S 5 was also integrable giving rise to the idea that 
the integrability structure of the string on AdS 5 x S 5 could be used as a vehicle 
to understand the AdS/CFT correspondence in detail, at least in the maximally 
super symmetric case ||12| . However, integrability is not a requisite of the AdS/CFT 
correspondence, so apart from circumstantial evidence for integrability at weak and 
strong coupling it is not clear that this will be the final answer to the AdS/CFT 
puzzle. 

Obviously, having an integrable structure is a strong constraint on the string 
dynamics and it has been shown that many interesting general perturbations of 
the maximally supersymmetric background destroy the integrability properties (see 
[ P~8| ] for example). In a similar vein, integrability can also happen for spin chain 
models with boundaries and it is interesting to examine if such integrable structures 
can also appear in the perturbative Yang Mills theory and use them as a tool to 
examine other configurations of matter in the AdS/CFT correspondence. Indeed, the 
natural boundaries for open strings are D-branes, so one can try to see if D-brane 
solutions preserve the integrability of the string theory and then try to solve the 
string spectrum to understand the D-branes, much like solving the string spectrum 
should help to understand gravity. 

Here, one can envision two different types of D-branes. The first corresponds 
to adding defects to the field theory so that there are "flavor branes" of various 
dimensionalities, and the quarks can act as boundaries for the spin-chain string. 



These branes are always infinite. This has been explored in |1J, [15], [16]], with the 
result that the one-loop set of anomalous dimensions are actually integrable. 

A second approach is to study finite volume D-branes on the AdS geometry. 
These correspond to non-perturbative states of finite energy in the string descrip- 
tion, so one can try to find the dual description of these states in the field theory. 
Because of issues of strong coupling physics, it is better if the corresponding states 
are supersymmetric to have some protection for the calculations. From this point 
of view, one has to look for supersymmetric D-brane configurations in the AdS ge- 
ometry, and to their corresponding dual description in the field theory. Once these 
D-branes are found, one should be able to describe the open strings attached to them 
and study if these strings can be described by a boundary spin chain model. After 
this is done, one can ask if the model corresponds to some form of integrability or 
not and if this integrability is of the familiar Bethe- Ansatz form or not. A lot more 
care is needed in this case because the branes are finite. This finiteness implies that 
the D-brane can back-react to the presence of the string and the dynamics can be 
much more complicated than in the case of an infinite brane. 
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This is the problem that we will concern ourselves in this paper. The branes 



under question are going to be giant gravitons [17|. These preserve half of the 



super symmetries and their dual field theory states are known very well [IIR O, 



How to add strings to giants was discussed in the works |21], g^, |23], |24j which also 
included a description of how the enhanced gauge symmetry of coinciding D-branes 
could be understood. One of the main difficulties in performing the calculations is 
that the dimension of the operator is of order N and it is harder to separate the planar 
and non-planar contributions. The dimension of the operator implies that there is a 
combinatorial enhancement of the usual 1/N suppressions, so it is not clear a priori 
that there is a well defined procedure that works in the general case, and one has 
to work example by example to understand the dynamics of the purported D-brane 
state. 

The one loop spectrum of anomalous dimensions for strings attached to a maxi- 
mal giant graviton was described in |25j , were it was found that the one loop planar 
anomalous dimensions correspond to an ordinary spin chain model with integrable 
Dirichlet-like boundary conditions. This work was extended to study what spin chain 



corresponds to a more general giant graviton in [26|, where we found that the spin 
chain in question has a variable number of sites and therefore it is not an ordinary 
spin chain model anymore. After a bosonization transformation, we found that the 
spin chain model could be also understood in terms of a system of a Cuntz oscillator 
chain model (a boson chain, where each spin corresponds to a single boson Fock 
space) with non-diagonal boundary conditions. The non-diagonal boundary condi- 
tions imply that the total boson number is not conserved. Our setup is also the 
limit q — > of a g-boson chain model with non-diagonal boundary conditions, where 
the q boson is defined by the algebra aa) — qa^a = 1, with q a real number. More 
recently, it has been claimed by p7| that in the maximal giant graviton case, there 
seems to be a problem with integrability beyond the two loop order, in contrast to 
the closed spin chain model |29|, |3(| . Agarwal argued that the Bethe Ansatz breaks 
down by direct calculation at two loops and therefore the model is not integrable. 
Another claim that studies the consistency of a Bethe Ansatz for the maximal giant 
graviton by Okamura and Yoshida suggest that the BMN limit breaks down instead 



P^ ] . We have a different interpretation of these facts: if the system is integrable, it 
will not realize integrability by a Bethe Ansatz. To do this, we study more general 
giant gravitons. 

In this paper we extend our analysis of |2(| to try to understand the full spectrum 
of the variable length spin chain model. The hamiltonian of the spin chain model is 
given by 

L L-l 



H = 2A 2_] ®i®'i ~ ^ /J(a[q;+i + a,ia] +1 ) 
1=1 i=i 

+2\a 2 + Xa(a\ + oi) + \a(a j L + a L ) , (1.1) 
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and the operators hi obey the Cuntz algebra for a single species, 

a 4 aJ = J, a\a l = /-|0)(0| , (1.2) 

where operators corresponding to different sites commute. l 

In our exploration of the model we will find very interesting phenomena that at 
first sight seem to contradict a very naive intuition about the system. The naive 
intuition is that the open string ending on the D-brane has a discrete spectrum, be- 
cause one has a finite energy configuration of matter in the AdS$ geometry. However, 
there is a technical issue with the fact that the energy is going to infinity as we take 
N large, because the mass of the D-brane scales with N. The naive answer, in spite 
of being technically correct, might also include some non-perturbative information 
on the finiteness of N, that is not necessarily reflected in the strict large N limit. 
Indeed, we find that the spectrum of the spin chain model is continuous. Upon 
thinking further about it, it becomes obvious that the spin chain answer correlates 
very well with the AdS geometric intuition and we will explain how this puzzle is 
solved in detail. 

The continuous string spectrum also has some other consequences for the integra- 
bility program. It seems to indicate that if integrability is present, it will not be of the 
form of a Bethe ansatz solution. From this point of view, if integrability is present, 



the conjectures in |30| need to be extended appropriately, or it is even possible that 
they might need to be revised completely. For finite size systems, a Bethe ansatz 
solution would predict that the spectrum of the system is related to solutions (roots) 
of polynomial equations of high degree, giving rise to a discrete spectrum of config- 
urations. Does this fact destroy the integrability of the open string? We don't think 
so. Although it is hard to come by with integrable systems that can not be solved 
by Bethe Ansatz, there are some well known examples of integrable systems that 
manage to be integrable without having a Bethe ansatz solution. Indeed, the c = 1 
matrix model is such an example. This model is dual to a two dimensional string in 
a linear dilaton background. It is very well known that the number of particle excita- 
tions is not conserved when scattering of the Liouville wall, so this non- conservation 
of the number of particle excitations is contrary to the typical Bethe ansatz solution. 
There are other hints that a Bethe ansatz is not the complete story for describing the 



individual string dynamics at strong 't Hooft coupling |3l| . We believe our results 
in this paper show that the Bethe ansatz breaks down perturbatively a lot faster for 
the open string than for the closed string. 

Our paper is organized as follows. We begin in section |2] by reviewing the physics 
of giant gravitons. We show how the local geometry near an open string spinning 



1 Since spin chain models also have applications to other areas of physics, the reader who is only 
interested in the analysis of the spin chain model can jump to section and should also read sections 
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around the moving giant graviton arises from a Penrose limit of AdS^ x S 5 . We 
quantize some of the bosonic modes of the open string in that background, which 
is valid for short strings that sense the local geometry. We then consider longer 
strings that carry two angular momenta, one of them in the same direction as the 
giant graviton. Then we expand the Polyakov action in the limit of large transverse 
angular momentum just like with closed strings J32], [M . 



In section || we move to the dual description in SYM theory. We review the 
derivation of the matrix of anomalous dimensions interpreted as the Hamiltonian of 
a spin chain with variable number of sites. We then map this problem to a lattice of 
bosons with sources and sinks at the boundaries. Then, we show numerical evidence 
for the agreement of the bosonic Hamiltonian spectrum with the plane wave spectrum 
of the dual open string. We also use coherent states to obtain a sigma model action 
in the limit of large R-charge. It agrees with the Polyakov action of the open string 
in the large momentum limit, in a particular gauge. Moreover we show how the SYM 
theory gives the correct boundary conditions for the open strings. 

In section |] we show that the spectrum of the anomalous dimension matrix con- 
tains continuum energy bands. We then argue that these can represent an instability 
of the dual D-brane where the open string can absorb a significant fraction of the 
angular momentum of the giant, making it collapse. In fact, we show that the lowest 
energy band is accessible to semiclassical string excitations. 

In section || we show some numerical evidence that lead us to conjecture the 
integrability of the Hamiltonian. Finally, in section |6|, we discuss our results and a 
list of open problems that hopefully will be of interest for the reader. Some of the 
details of the calculations for the different sections are shown in the appendices. 

2. Open Strings on giant gravitons in AdS§ x S 5 

In this section we study giant gravitons in S 5 and their open string excitations from 
the point of view of string theory /supergravity. In the next section we will see how 
this stringy description arises from the dual SYM theory. 

We begin by briefly reviewing some of the basic geometrical properties of the 
spherical D3-branes existing in the background of AdS^ x S 5 that wrap an S 3 inside 
the S 5 ||17|| . These spherical branes are commonly called giant gravitons. 

They are 1/2 BPS and their classical stability is due to the presence of the five- 
form flux which exactly cancels their tension. As we will see, the movement of the 
non-maximal giant gravitons gives rise to non-trivial boundary conditions to open 
string excitations. We are interested in studying fast moving (nearly supersymmetric) 
strings on the D-brane (close to the speed of light and near the point particle limit), 
as for these strings the dynamics simplifies. These should be closely related to null 
trajectories in AdS$ x S 5 that are contained in the D-brane worldvolume trajectory. 
Moreover, as we now show, for the special supersymmetric case the geometry near 



-6- 



this null trajectory is the usual pp-wave of type IIB supergravity, so we will focus on 
this case. 

2.1 Spherical D3-branes in AdS 5 x S 5 

Spherical stable branes are known to exist in the maximally supersymmetric back- 
ground of type IIB supergravity AdS$ x S 5 [O. Indeed, these branes are actually 



supersymmetric. They can wrap an S 3 either inside the S 5 or inside of AdS$ |J5|. J36fl . 
In this section and in the rest of the paper we will focus on the former ones. 
To begin with, let us write the metric of the AdS§ x S 5 background as 

ds 2 = R 2 {- cosh 2 pdt 2 + dp 2 + sinh 2 pdVt'l + d9 2 + cos 2 9dip 2 + sin 2 9d^\) , (2.1) 

and the 3-sphere metrics as 

cZO| = dip 2 + cos 2 ipdr] 2 + sin 2 pd£ 2 , 

dVt'l = dip' 2 + cos 2 p'dT]' 2 + sin 2 p'di' 2 . (2.2) 

In these coordinates, the Ramond-Ramond 4-form potential takes the form 

C (4) = 4vr iVc/ 2 (sinh 4 pdtAQ' 3 - sin 4 9 di\) A fi 3 ) , (2.3) 

and we have the relation R 4 = Airg s Na' 2 . 

McGreevy, Sussking and Toumbas found spherical D3-branes carrying angular 
momentum in the S 5 ||17|| . They are supersymmetric solutions of the brane action, 
that expand in the 3-sphere f^. Choosing the static gauge, the parametric coordi- 
nates of the brane (<t , cr 1 , a 2 , er 3 ) can be identified with space-time coordinates 

t = <T Q , (p = 0"! , X] = 02 , i = CT 3 • (2-4) 

The brane carries angular momentum along the ip direction and it is located at p = 
and at a constant 6 = 0q. 

^ = #7 ), p = o, e = e , (2.5) 

The equations of motion are solved with ip constant. More precisely, and indepen- 
dently of 6q, by ip = 1. Thus, the center of mass of the giant graviton is moving 
along an equatorial null trajectory. However, each element of the giant is moving in 
a time- like orbit of radius r e \ = Rcos9 . The radius of the giant also depends on 9 , 
r gg = Rsin9 . 

For ip — 1, the momentum conjugate to ip becomes 

p = Nsm 2 9 . (2.6) 
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We can use this relation for highlighting a relevant feature of giant gravitons: their 
radii grows as they increase their angular momenta 



r 99 = R \jjj- ( 2 - 7 ) 

Since we are considering giant gravitons that expand in S 5 , the radius is bounded by 
r gg < R and hence the angular momentum is also bounded by the number of units 
of five-form flux on the sphere 

p<N. (2.8) 

When the equality is satisfied the brane solution is known as maximal giant gravi- 
ton. Notice that in this case the radius of the orbit of an element of the giant 
r e i = a/1 — p/N shrinks to zero and all the angular momentum comes from the 
term with the Ramond-Ramond form in the brane action. Being the maximal giant 
graviton static is the reason why this special case turns out to be simpler. However, 
considering non-maximal giant gravitons gives rise to interesting and novel phenom- 
ena. 

2.2 A Penrose limit for an open string on a non-maximal giant 

We obtain in this section, the effective geometry seen, in the large N limit, by an 
open string with angular momentum J growing as y/N, attached to a giant graviton 
whose angular momentum p is growing proportionally to N. Then, the angle 8q is 
kept constant and the radius of the giant diverges. The open string will be effectively 
attached to a flat D3-brane in a pp-wave background. For the maximal giant graviton 
this was already analyzed in J22|. We now focus on smaller or non-maximal giant 



gravitons. 

As we already said, each element of the giant graviton is moving in a time-like 
trajectory. In order to travel in null trajectory, an observer on the giant should 
be spinning fast along it. Then, we will study open strings with two components 
of angular momenta. In particular, we now consider a trajectory along ip and r), 
keeping ip = 0. A null trajectory should satisfy 

R 2 (-i 2 + cos 2 6 Q i) 2 + sin 2 6 f] 2 ) = . (2.9) 

Since t — 1 and ip = I, one necessarily has fj = ±1. Although it is not obvious at 
first sight, this null trajectory is a null geodesic of A0IS5 x S 5 and then a customary 
Penrose limit can be taken. At this point it is worth highlighting the ratio of the 
angular momentum in both angular directions 

4=„t»* = £(.-£). (2,0) 

We will later see the natural emergence of this quantity in the gauge theory side. 



By following the standard Penrose limit procedure (see appendix [A] for details) 
we get a plane-wave geometry with the metric 

6 6 

ds 2 = —Adudv + Aydudx — 2_. z^du 2 + dx 2 + dy 2 + 2_. dz 2 . (2-11) 

a=l a=l 

The RR 5-form field strength becomes 

F( 5 ) = 2du A (dzi A dz 2 A dz 3 A dz 4 + dz 5 A dz% A dz 7 A dz s ) . (2.12) 

This is nothing but the usual maximally supersymmetric pp-wave of type IIB su- 
pergravity || displayed in unusual coordinates. This can be explicitly seen with an 
appropriate coordinate transformation |J7], |58 |. 



2.3 Open Strings on the pp-wave Geometry (Short Strings) 

Let us now consider the open string theory in the pp-wave geometry Q2.11| ) corre- 



sponding to an open string spinning along a non-maximal giant. This is a good 
description for short strings which sense the local geometry near the giant graviton. 
We will focus on the bosonic sector of the superstring and, in particular, on those 
modes that later will be compared to the gauge theory predictions. More precisely, 
we will focus on open strings carrying two angular momenta on the sphere with one 
of them in the direction of movement of the giant graviton. 
The bosonic part of the superstring action is 

S = / dr / da (—Ad a ud a v + Ayd a ud a x — z i z i d a ud ol u 

4vra' J J 

+d a xd a x + d a yd a y + d a z l d a z i ) . (2.13) 

Fixing the light-cone gauge with the usual procedure 

u = 2a'p u r, (2.14) 

we obtain the light-cone action 

Sic = dr da (d a xd a x + d a yd a y — Amyx + d a Zid a Zi + m 2 ZiZi) , (2.15) 

4vra' J Jo 

where dots and primes refer to derivatives with respect to r and a respectively. We 
have also defined the mass m = 2a'p u . As in the case of closed strings [[39|j , this action 
can also be derived as the quadratic order expansion of the string action around the 
classical solution 6 = 9 and if) = t. 

The equations of motion from (|2.15|) are 



{d a d a - m 2 ) Zi = , (2.16) 

d a d a x + 2my = , (2.17) 

d a d a y - 2mx = . (2.18) 
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For Zi we obtain the usual massive equations. For x and y it is convenient to define 
a complex field w = x + iy which satisfies 

d a d a w - 2imw = , (2.19) 

A solution of ( J2.19 ) is w = e~ tmT W whenever W is a solution of the massive equation 

(d a d a - m 2 ) W = . (2.20) 

We now have to specify the boundary conditions for these fields. To this end, we have 
to keep track of the D-brane position in the Penrose limit. This can be done using 
the coordinate transformation of the Penrose limit (see appendix [B]). The original 
boundary conditions are translated into Neumann boundary conditions for u, v, Z\, Z2 
and Dirichlet boundary conditions for x, y, £3, £4, z 5 , z 6 . 

The mode expansion and the canonical quantization of fields z±,...,Zq (those 
satisfying massive equations of motion), goes exactly as in [10 1. We concentrate on 



excitations of fields x and y which later will be related to the dual description. 

Solutions of (|2.20|) satisfying Dirichlet boundary conditions can be expanded as 



W = J2 Mncr) J— (&e-*"* T + foe*""*) , (2.21) 

n>0 ' UJn 

where u n = +y/n 2 + m?. The factor ^jAa'/uj n in the coefficients is included for later 
convenience. The expansions for the original fields x and y are 



x 



~ J2 Mno-) J^f (P n e-^ T + (3 n e- iujir + P* n e^ T + (3* n e^^ , (2.22) 



n>0 ' "" 

i sr-^ . , , /4a 
2 



^sin(na)W— (/3 n e-^ r -/3 n e-^ r -/3:e^ T + /?:e^ T ) , (2.23) 

n>0 ' n 

where now, u^ = u n ± m. The string spectrum can be obtained by canonical 
quantization. With our normalization, the coefficients upgraded to operators satisfy 
two set of mutually commuting oscillator-like algebra, 

\p n , Pi] = 6 nm , \J3 n , (3 m \ = \pl Pi] = , (2.24) 

n , $l\ = 5 nm , n , p m ] = 01 fa] = . (2.25) 

The spectrum of the sector we are considering is obtained by acting with the creation 
operators /3^ and p^ on a vacuum state satisfying P m \0) = $ m \0) = 0. The light-cone 
Hamiltonian can be expressed in terms of the number operators corresponding to the 
oscillator operators. Scaling r and o by 2a'p u 

"2na'p u 



1 pznap- 

HZ-j^l d a(x* + f + X " + V « 

l — J2 {<Mn + W+APn) ■ (2.26) 



2a'p 

* n>0 
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We see that excitations created by /?t have more energy than those created by p^. 
Expanding the square roots in u n 



n>0 L * \ f / 



(2.27) 



where iV n = ^/3 n and iV n = ^/3 n . 



The Hamiltonian and angular momentum generators are, in the light-cone gauge, 

H lc = - Pu = i—, (2.28) 

P"=-^ = ~- (2-29) 



Using the change of coordinates ( |A.1| ) we can express them in terms of the original 
generators 

i ( d d d \ A + Jjf, + Jrj 



P 2R 2 \dt di) dr]J 2R 2 ' ^' 3 ^ 

Looking forward to a gauge theory interpretation we can label the angular mo- 
menta as (c.f. ( |2.10D ), 

J V = L, (2.32) 

a 2 

J^, = cot 2 9 L = -L . (2.33) 

1 — cr 



where we define a = a/1 — p/N. Then, the sum of angular momenta is 

J v + J^ = —^ ■ (2.34) 

I — cr 

For a finite light-cone energy it is required that A ~ L/(l — a 2 ) and then p u ~ 
L/R 2 (l — a 2 ). Finally, using R 4, = Ang s Na' 2 and A = g s N/2n, the energy of each 
excitation is 

Xn 2 (l-a 2 ) 2 n 2 \n 2 (l - a 2 ) 2 n 2 

E n « ^ [i L2 a)n , E n ^2+ An[i L2 a)n , (2.35) 

for x and y respectively. Thus, as expected, we see that string theory predicts a 
BMN limit for the anomalous dimension of the dual operators describing these open 
string excitations. We will come back to the dual interpretation of these energies in 
section 3. 
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2.4 Semiclassical Open Strings (Long Strings) 

If we want to consider more general open strings ending on the giant graviton, we 
need to include the full AdS$ x S 5 background in the Polyakov action. However, it is 
well known by now that we should expect a classical description in the large angular 
momentum limit 2 . Again, as we did for short strings in the pp-wave geometry, we 
focus on the excitation of two coordinates of the string. More precisely, the two 
coordinates of the sphere 5* 5 subject to Dirichlet boundary conditions. The reason 
is that the description of these excitations in the dual theory is more easily isolated 
from the rest. Moreover we will fix a particular gauge of reparametrization invariance 



different from the standard conformal gauge, following closely |34j . This election will 



be the reflection of the a particular labeling of the operators in the dual gauge theory. 



These results were already presented in a previous letter |2(| . We reproduce them 
here providing more details of their derivation. The starting point is the Polyakov 
action in phase space. We can write the conjugate momenta as, 

Pfl = -G^Aaox" + Bdix») , (2.36) 

where A = y/—gg 00 , B = y/—gg Q1 and g ab is the worldsheet metric. The Polyakov 
action then takes the form 



where, 



S P = v^ J drJ*^C, (2.37) 



C = p„d x^ + X -A' x \G^ VixVv + G^d^dxx"] + BA-%0^ 1 . (2.38) 



Here we have factorized the radius of AdS§ and S 5 so that by the AdS/CFT cor- 
respondence X Y m = gYM-N = 8n 2 \ = R 4 /a 12 . Moreover, A, B play the role of 
Lagrange multipliers implementing the constraints T a b = 0. 

For an open string traveling with the giant graviton and with excitations on the 
sphere only, the effective geometry is M x S 5 . We can write the corresponding metric 
as. 

ds 2 = -dt 2 + \dX\ 2 + \dY\ 2 + \dZ\ 2 , (2.39) 

where |X| 2 + |y| 2 + |Z| 2 = l. The giant graviton will be orbiting in the Z direction 
with Z = a/1 — p/Ne lt and will wrap the remaining S 3 . We will put our string at 
X = 0. Then we define the coordinates, 



'■ 



Z = re %(t -* } , (2.40) 

Y = ±Vl - r 2 e iip , (2.41) 



The literature on this subject is very extensive, but for a nice review see 42 



- 12- 



for which the giant graviton is static at r = yl — p/N and, say, = (we can 
always shift by a constant). The metric becomes, 

ds 2 = -(1 - r 2 )dt 2 + 2r 2 dtd(p + - -dr 2 + r 2 d(p 2 + (1 - r 2 )dip 2 . (2.42) 

The connection between these coordinates and the plane wave coordinates is shown 
in appendix |A|. 

The momentum in ip is conserved and is given by, 

PIT 1 

L = v 7 ^ / ttP<p = V^ymJ ■ ( 2 - 43 ) 

Jo 21T 

We choose a gauge that distributes the angular momentum p v homogeneously along 
the string. Furthermore, we choose r to coincide with the global time in the metric. 
Thus, our gauge is 

t = T , Pip = 2J = const. (2.44) 



These type of gauges were introduced in [ 33 1 



We then implement the constraints that follow from varying A and B in ( |2.37| ) 
directly in the action. This allows us to write the Lagrangian in terms of the momenta 
p r and ptf, and the fields r and and their derivatives. These are the two fields subject 
to Dirichlet boundary conditions. We get (up to a total derivative in r), 



C = (1 - (fyp^ + rpr - J apl + bp 2 r + 2cp r p 4> + d , (2.45) 
where, 

«= 1 + ( 1 -0(^ + ^) , (2.46) 

^ = (l-0(^ + l) , (2-47) 



r>'(l - r 2 ) 

AJ 2 
4 ^ 2 ■ 2,2 , r 



(2.48) 



,2 



d=-^— 2 +r 2 <P' z + - -. (2.49) 

1 — r l 1 — r z 

As usual, dots and primes denote derivatives with respect to r and a respectively. 
The remaining gauge freedom can be fixed by demanding that the equations of 
motion for p r , p^, r and <fi that follow from ( |2.45| ) agree with the ones derived from 
the original action ( 2.37|) p4[ . This will fix the value of B in terms of p r , p$, r and 

0- 
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Since the momenta p r and p^ enter the Lagrangian ( 2.45 ) algebraically, we can 
solve for them using their equations of motion and write the action in terms of the 
fields r and <f> and their derivatives. We get, 



£ = - 



d 



1/2 



1 - a v\ - bPr + 2cp <l) p r 



1 + (1- 



— rp r 



(2.50) 



where, 



re — (1 



Pr 



ab — c 2 
fa — (1 — d>)c 



(2.51) 
(2.52) 



ab — c 2 

Finally, we take the limit J — > oo and assume that the time derivatives are of 
the order Oqx^ ~ 1/J7" 2 . This last condition can be made more precise by solving 



for the time derivatives in terms of the spatial derivatives as in | 34fl . However, in 
our case this is not necessary because we just want the lowest order in 1/ J . Thus, 
rescaling o — > na, we get to lowest order, 



S 



-L 



i 

dt J da 

o 



r 2 <j) 



+ 



1? 



r' 2 + r 2 6' 2 ) + O 



L 4 



(2.53) 



We note that L serves as an inverse "Planck constant" and so L — > oo corresponds 
to a classical limit as promised. This coordinates are subject to Dirichlet boundary 
condition, which are expressed as, 



k=0,l 



|cr=0,l 



N 



const. 



(2.54) 
(2.55) 



3. Open Strings on Giant Gravitons from J\f = 4 SYM 

In this section we study the dual description of open strings on giant gravitons given 
by M = 4 SYM. It is very interesting to see how the complete geometrical picture 
discussed in the last section is recovered from the gauge theory. In particular, we show 
the emergence of the BMN spectrum, the Polyakov action, the Dirichlet boundary 
conditions and the reparametrization invariance of the string world sheet. All of 
these geometrical ingredients are encoded in the matrix of anomalous dimension of 
the operators dual to giant gravitons with open string excitations. 

There is considerable evidence that the dual operators to these D3-branes are of 
determinant-like form 0, |19|, |23| (see appendix for our conventions) 






/r\p __ Jlf'Jp y%l _ _ _ y l v 



.11 



■IP 



(3.1) 
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where Z is one of the three complex scalars of M = 4 SYM. These are chiral operators 
and, by supersymmetry, their dimensions are determined in terms of their R-charge: 
A = J (BPS condition). For O p we see that A = p < N, and hence this operator 
obeys the "momentum" bound of the giant gravitons. This makes sense from the 
AdS/CFT correspondence since we should identify the U(l) charge of ( |3.1| ) as the 
angular momentum of the D-brane and A as its energy. 

The operator dual to a single string attached to an S 5 giant graviton of momen- 
tum p is obtained by appending a word to the determinant-like operator |E| 24] , 
which represents the open string 3 , 

where W is a "word" built out of the scalar fields X, X,Y,Y, Z, Z, as long as we 
only consider excitations on the S 5 directions. Each scalar field carries a unit of U(l) 
charge corresponding to a unit of angular momentum in one of the three orthogonal 
planes that cut S 5 in the dual string picture. 

We also have the condition that the Z field is not allowed to be at the borders of 
W [24"]. This constraint comes from the fact that when a Z is at the border of W, the 



operator factorizes as a bigger giant plus a giant of the same size with a closed string 
(see appendix |C]). In fact, in general, the operator ( |3.2fl can be expanded in terms of 
traces ||23| . However, for p sufficiently large the mixing of ( |3.2| ) with closed strings 
(traces) is suppressed in the large N limit. For this reason we will only consider 
operators with yN <p<N. 



3.1 Open Strings as Variable Length Spin Chains 

We are interested in computing the mixing matrix of anomalous dimension for this 
type of operators at the one-loop approximation and large N limit. Thus, we will only 
consider planar diagrams. As usual, we begin by defining the correlation function 



M AB = (O A (x)O B (0)) , (3.3) 

\ / free + interacting 

where the operators have been normalized according to 

6 A {x) = ° a{x) 1/2 , (3.4) 

and where A, B are collective indices labeling different giant graviton configurations 
with a single string. At one- loop and in the large N limit, M AB will have the following 
general form, 

M AB = r-^(5 AB - 2 T AB log(|x|A) + ...), (3.5) 



3 More precisely, as we will see later, the dual operator to the classical non-maximal giant graviton 



is actually a coherent state of operators like (3.2) with different values of p 
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where Ao is the classical dimension, Tab is the matrix of anomalous dimension and 
A is an ultraviolet cutoff. We then identify the anomalous dimension matrix with 
the Hamiltonian of the corresponding string quantum states, 

T AB = (iPa\H\iPb) , (3.6) 

where Oa — \iPa)- This is just the operator state correspondence which is available 
for any CFT: the Hamiltonian corresponds to taking radial time and compactifying 
the CFT on a round sphere. 

In |25 1 the case of a maximal giant graviton was considered and it was shown 



that the resulting planar anomalous dimension matrix corresponds to an integrable 
open spin chain with SO (6) symmetry. For the case of non- maximal giants, most of 
the field contractions are the same as with the maximal giants, the difference being 
the particular combinatorics of the e symbol for p ^ N. However as we remarked 
in ||26|| , there is a very important new interaction in the case of the non-maximal 
giant graviton: Z fields can be exchanged between the word W and the rest of the 
operator. This accounts for the exchange of angular momentum between the giant 
and the string since the open string is "dragged" by the non-trivial movement of the 
non-maximal giant graviton. 

The details of the combinatorics involved for this and the rest of the interactions 
between the fields in the operator are explained in the appendix |C[ Here we just 
quote the main result focusing on a SU{2) subsector involving the fields (say) Y 
and Z 4 . There are two kind of terms in the mixing matrix of anomalous dimension. 
Those corresponding to bulk interactions in the spin chain Hamiltonian, which where 
already present in the maximal case. They are nothing but the bulk terms of a 
SU(2) XXX spin chain. The non-maximality of the giant gives rise to new terms 
proportional to the quantity a = a/1 — p/N. The most important is a boundary term 
for the spin chain Hamiltonian, coming from the correlation function representing 
the exchange of a Z field between the word W and the determinant 



irfp+u, n fyg / n x X 2Alog(A|x|) I p W = Yw, 

(°w (*)<^(°)> ^— V A 7 ' Wlth ' W' = YZw. (3 ' 7) 

This correlation function introduces a variability in the spin chain length at the 
boundaries. However, it is not a priori clear how to deal with this SU(2) spin chain 
of variable length. In the following we present how, by relabeling the operators, we 
can translate the spin chain of variable length into a more manageable bosonic lattice 
of fixed length and with sources/sinks of bosons at the boundaries. 

The most general word in this sector is completely specified by the number of Z 
fields between two consecutive Y fields, 

(YZ ni YZ n2 Y . . . YZ nL Y){ = \n h n 2 , . . . , n L ) . (3.8) 



4 The labeling of the scalar fields by Y and Z is done on purpose to make the connection with 
the coordinates used in section 2.3 more explicit. 
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Thus, the word has a fixed number L + 1 of Y fields, representing L + 1 units 
of angular momentum in Y direction in the string. The number of Z fields and the 
corresponding angular momentum is measured by the number operator n = ^TZj. 
The size of the giant graviton is measured by p = (p + 1)1 — n, where p is the (fixed) 
total number of Z fields in the operator. The exchange of Z fields occurs only at the 
first and last site of the bosonic lattice. We will see that the choice of leaving fixed 
the number of Y fields in our labeling is in correspondence with the gauge choice 
fl2.44|) in the Polyakov action. 



Since there is equal probability of a Z entering or leaving the word (see appendix 
|CJ) , one expects that for sufficiently large p the lowest energy states have ( {p) —p) /N ~ 
in the large N limit. That is, the backreaction to the size of the giant graviton 
should be negligible. Therefore, in what follows we assume that p ~ 7iV with 
< 7 < 1, and approximate p/N m p/N in all matrix elements of the Hamiltonian. 
Later, we will see that this is a consistent approximation. 

The Hamiltonian dual to the anomalous dimension matrix for this sector then 
takes the form 

L L-l 



H = 2A y^ ajai — X y^(alai + i + d/d]^ 
i=i i=i 

+2A a 2 + X a (a\ + di) + A a (a L + a^) , (3.9) 



where L < yN and the operators hi obey the Cuntz algebra for a single species [El] . 



aia\ = I, a\a t = I -\0)(0\ , (3.10) 

where operators corresponding to different sites commute. 

For our purposes it will be more useful to think of this algebra as the q — > limit 
of the deformed Weyl algebra, dd^ — gd^d = 1, [h, d^] = a\ [n, a] = —a. The number 



operator ft can be constructed in terms of the oscillator operators as in [[44 . 

Terms in the first line of the Hamiltonian (|3.9|) come from the bulk interactions 
that were already present in the maximal case. The first term indicates that each 
bosonic oscillator contributes with an energy 2A whenever its site is occupied. The 
second term is a hopping interaction for bosons to move between sites, so that the 
energy is reduced with bosons which are not localized. The second line of ( |3.9| ), 
apart from the constant term, provides source and sink terms at the boundaries of 
the bosonic lattice, which give rise to non-diagonal boundary conditions, since the 
total boson occupation number does not commute with the Hamiltonian. 

In the next subsections we will discuss two different sectors of the Hamiltonian 
( |3.9| ) and how they give rise to the BMN limit of short strings and the semiclassical 
limit of long strings discussed above. 



-17- 



3.2 Evidence for a BMN Limit 

At this point, it could be a significant verification of the validity of our dual descrip- 
tion to see the BMN spectrum ( |2.35| ) arising from the Hamiltonian ( |3.9| ). Unfortu- 



nately, we have not been able to diagonalize the Hamiltonian ( |3.9| ). Nevertheless, 
evidence for its integrability will be given in section 5. 
However, we do know the ground state: 



CO 



|*o> = (1 - « 2 ) L/2 E (-«r + "' +nL \ni,...,n L ), (3.11) 

ni,...,ri£=0 

and it has energy E = 0. The expectation value of the number operator for the 
ground state is, 

<*o|A|*o> = ^ (l - £) , (3-12) 

which is generically of order L, unless p << N. Since L < \fN we see that the 
backreaction to the giant is indeed small compared to p ~ N. Moreover, note that 
setting a = gives the familiar ferromagnetic ground state of the maximal giant 
graviton: 

lim |tf o> = |0, 0, . . . , 0) ~ (YY ■ ■ ■ Y){ . (3.13) 

The expectation value ( J3.12 ) gives the amount of angular momentum that the 
string acquires in the direction of the movement of giant. Dividing by the fixed 
number of Y fields we obtain the ratio of angular momentum in the two directions 
of the S 5 . 

Jf ^ number of Zs in W ^ (^o\n\^/ ) N (^ p\ , . 



Jrj number of Ys in W L p V iV 

This is precisely the same ratio of angular momentum components of the null-geodesic 
we used to take the Penrose limit ( p.lOj ). Thus, as expected, the ground state ( 3.11 ) 



corresponds to a point like string traveling with the giant, and small fluctuations 
around it should correspond to the modes of the open string in the pp-wave back- 
ground. Unfortunately we cannot solve for the lowest energy modes of (|3.9|) in gen- 



eral. However, we can treat perturbatively the boundary terms (a = y/l — p/N < 1) 
and try to reproduce the spectrum ( ^.35| ) at lowest order in a. 

Before doing this, a few comments are in order. First, let us look at the case 
of the maximal giant graviton (a = 0). In this case there is a precise dictionary 
between the plane wave excitations and the corresponding gauge theory operators 
] . The ground state ( |3.13| ) is excited by adding Z fields to the word W with some 
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momentum determined by the boundary conditions. For example, the first excitation 
is given by a single Z: 



m = yf^ X>in (0i) 10 > where 10 - (Y l ZY L+l - l ){ . (3.15) 



i=i 
The anomalous dimension of this operator is, 



E^ = 2A 



1 — cos 



L + l 



-^-, (3-16) 



On the other hand, the first excited state of the open string in the pp-wave 
background has energy, 

\ii 2 n 2 

^n(String) « 1 + — j^~ ■ (3.17) 

The factor of 1 can be interpreted from the increase of the classical dimension of 
the word by the insertion of a single Z field. Then, the anomalous dimension ( |3.1G| ) 
agrees with the small correction in fl3.17| ). 

However, for the case of the non-maximal giant graviton there is no precise 
notion of the length of the word, and hence its classical dimension. This is because 
Z letters can enter and leave the word. This means that we do not have a clear cut 
distinction between what we call the "string" and the "giant". In fact, we only know 
the average length of the word to leading order in L: L + 1 + (n) ~ O(L). Adding 
or subtracting a finite number of Z fields to our operators will have no effect in the 
anomalous dimension in the planar approximation. 

In the spectrum of open strings on the pp-wave geometry we found what appears 
to be two independent energy modes (|2.35|) differing by a factor of even integers. 
But this is precisely the ambiguity in the classical dimension of the corresponding 
operators. Therefore, we should only compare the anomalous dimensions of our 
operators to the 0(X/L 2 ) corrections in the string theory spectrum. In fact, we see 
that setting a = in (|2.35|) gives precisely the anomalous dimension (|3.16|) . 

We now want to calculate the lowest order correction in a to the energy of the 
lowest BMN excitations of the maximal giant. For the sake of simplicity, we focus 
on the eigenstates with one-boson; that is, a single Z on the word. We write the 
Hamiltonian as 

H = 2\a 2 + H + aV , (3.18) 



where 



and 



L-l 



H = 2A ^2 a\h ~ A S ^2(a\a l+l + aih\ +l ) , (3.19) 

i=i i=i 

V = \{a\ + dli + 4 + ®l) ■ (3.20) 
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The ground state \ip°) of the unperturbed Hamiltonian Hq has zero occupation num- 
ber and zero energy (c.f. ( p,13|) ). The one-boson eigenstates and eigenenergies were 
given in Eqs. ( |3.15| ) and ( |3.16| ). 



The corrected energy is 



E„ 



2Xa 2 + Ei 0) + aV nn + a 2 



E 

kj^n 



IK 



nk 



E^ - E™ 



0{a 3 



(3.21) 



Matrix elements of V in the unperturbed basis are non-vanishing only for eigenstates 
that differ in one boson. Then, V nn is zero and the first correction is order a 2 . To 
compute the correction to the one-boson eigenvalue at this perturbative order, the 
only matrix elements needed are 



Mv\^°) 



A 



L + l 



sin 



nir 
L + l 



and MV\^ 2 J , (3.22) 



where |^) are the unperturbed two-boson eigenstates (in this case m is a collective 
index). However, it is hard to write down the two-boson eigenstates in closed form. 
One can not use a dilute gas approximation, because this would only be valid for 
small energies and we need \ipm) f° r an energies. Still, it is possible diagonalize 
numerically Hq in the two-boson subspace, i.e. to solve the eigenvalue problem 
for a L(L + l)/2 x L(L + l)/2 matrix, for different finite values of L and use the 
results to compute the energy correction. Then, we can see that as the number of 
sites increases, the correction not only captures the 1/L 2 dependence but also the 
proportionality factor tends to the one predicted by the plane-wave spectrum 



E ~ 

J - J n — 



A 



n 2 n 2 



2a 2 ). 



(3.23) 



For instance, in table 1 we show the shift of the n = 1 one-boson eigenvalue divided 
by the predicted shift (|3.23|) . 



Table 1: Shift in n = 1 one-boson energy. 



L 


ABi 


2tt 2 /(L+1) 2 


15 
20 
25 
30 


1.122 
1.091 
1.079 
1.060 



If our Hamiltonian is integrable (see section 5), then maybe some generalization 
of the Bethe Ansatz can be used to diagonalize it and provide the final proof of 
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the presence of the BMN limit. Integrability is not strictly necessary. There are 
other forms of obtaining the BMN limit exactly at strong coupling using different 
techniques ]45, 



>T], [PJ, see also J4?|. We also have uncovered some other qualitative 
features of the spectrum including the presence of continuum bands by doing pertur- 
bative diagrams that cast doubt on a solutions of the problem in terms of a Bethe 
Ansatz. We discuss these in section 13. 



3.3 The Semiclassical Limit 

In general, it is possible to obtain a semiclassical sigma model action governing 
the dynamics of a given spin chain. In |!2] it was shown that the sigma model 
obtained from the spin chain associated to the SU(2) sector of SYM coincides with 
the Polyakov action describing the propagation of closed strings in AdS^ x S 5 , in a 
particular limit. This comparison was extended to other sectors j|j|, ^9[ and to other 
realizations of the AdS/CFT correspondence [|5i 



1M- 



Now we will obtain a semiclassical sigma model action that governs the dynamics 
of bosonic lattice in the L — ► oo limit, using the coherent states basis for the path 
integral representation of the evolution operator. This was already discussed in our 
previous note [F2BJ. For completeness, we will use the g-deformed Weyl algebra and 
will set q — > when needed. In appendix fFJ, we give the definition of the coherent 
states for this algebra. 

One can then put a coherent state at each site of the bosonic lattice and construct 
the propagator between the coherent states in the usual way (see |k|). The resulting 
action is, 



S- 



/ dt li(zi . . . z L \ — \z 1 . ..z L )- (zi . . . z L \H\z 1 . . . z L ) j 

L 

^2fg( r i)<fii -2A (a 2 + a(rxcos0i + r L cos0 L )) 

L 

-2A ( Y^ r i - Yl TiTi+l cos( ^ ~ ^ +i )) 



dt 



,i=l i=l 



(3.24) 



where we defined, 



fq(x) 



x 



2 exp q (x z 
exp Jx 2 



(3.25) 



In the large L limit, the complex modulus r and the complex argument <fi become 
functions of a continuous variable < a < 1. We consider the action (|3.24|) in the 
limit L — ► oo and A — ► oo but keeping X/L 2 fixed and small. The result is the 



- 21 - 



following sigma model, 

S = -L 

-A 

-A 



1 

dt I da 

o 



/ 9 (r)0+^(r /2 H-rV 2 ) 



dt[ 
dt[ 



( 2 sin 2 d> + (a cos 6 + r ) 



2 • 2 

a sin - 



+ (acos0 + r) ] 



(T=0 



a=\ 



(3.26) 



where dots and primes refer to derivatives with respect to t and a respectively. If 
we take the limit q — > 0, the function above reduces to f q (r) — > r 2 /(l — r 2 ) and we 
see that the bulk action of ( |3.24| ) coincides with the Polyakov action in the large 



momentum limit ( 2.53 ). This gives a direct geometrical meaning to the fields r and 
of the coherent states: as spacetime coordinates. 
The classical Hamiltonian of the coherent states is 
•l 

'2 



(H) = t 



da{r' 



„2i/2\ 



+ x[ 



( 2 sin 2 4> + (a cos + r) ] | ^ =Q 



+ A [a 2 sin 2 (ft + (a cos <p + r ) ] | =1 



(3.27) 



One sees that the boundary terms will give rise to a large anomalous dimension of 
order ~ A unless, 



H<r=0,l 

0U=O,1 



a = 

IT . 



1 



P_ 

N ' 



(3.28) 
(3.29) 



These are exactly the Dirichlet boundary conditions ( |2.54j ) on open strings on non- 
maximal giant gravitons! The fact that the boundary terms give rise to a large 
anomalous dimension is nothing but the statement that moving the D-brane takes a 
lot of energy. 

In the limit L — > oo we get the following classical equations of motion: 






A 










. 



(3.30) 
(3.31) 



Looking at the Hamiltonian ( |3.27|) we see that the classical solutions to these equa- 
tions will have energies ~ A/L, which represent a small multiplicative correction to 
their "bare" energies: L + 1 + (n) + 0(\/L) ~ L(const. + 0(X/L 2 )). This is already 
a familiar result for closed semiclassical strings ||42|| . 
The average number of bosons in the lattice is, 



(n) 



1 r 2 



L / da- 



(3.32) 
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This provides a way of measuring length of the spin chain in the original XXX model 
formulation. In general, (n) is not conserved and therefore the string will oscillate 
in length. Using Eq. ( |3.30| ), this variation can be put simply as, 

d t (n) = 2± (l - |p) (0^=0 - 0W) • (3.33) 

Note however that we must ensure that (n) is bounded so that neglecting the back- 
reaction to the giant remains a good approximation. 

Moreover, we see how the particular choice of spacetime coordinates and world 
sheet gauge are encoded in the SYM side. We labeled our states using the word 
W only, which is dual to the open string. This is translated to the string side by 
choosing a coordinate system for which the giant graviton is static. Finally, by 
labeling the states as in ( |3.8|) we are explicitly distributing the angular momentum 
in Y uniformly along the string. This has very strong implications in the AdS/CFT 
correspondence, because we are seeing explicitly the reparametrization invariance of 
the string worldsheet: the gauge that makes the calculation more natural is different 
than the one considered in other semiclassical setups [^, Q. 



4. A D-brane Instability? 

He have seen that the anomalous dimensions of the dual operators of non-maximal 
giant gravitons with open strings excitations is described in terms of a spin chain of 
variable length. We offered an alternative description in terms of a bosonic lattice. 
There, the variability of the spin chain's length was translated into the variability 
of the boson occupation number. In this section, we study the possibility of having 
configurations with the occupation number growing monotonically in time. More 
precisely, we show evidence for the presence of continuum energy bands in the spec- 
trum of the Hamiltonian ( |3.9|) . If a configuration has enough energy to access the 
continuum, a more general time dependence than an oscillation is allowed for its 
mean occupation number. 

In particular, we argue that in the large L limit there is at least one band 
that is accessible to classical long strings with energies ~ X/L. This has profound 
consequences for the stability of the D-branes because the word W in (|3.2|) can be 
excited in such a way to absorb a large number of Z fields from the giant graviton. 
In the dual picture, this means that the string becomes longer and longer absorbing 
more and more angular momentum from the D-brane. A long string attached to the 
moving giant graviton, could suffer centrifugal forces, consequence of a non-geodesic 
movement, and if the string is long enough its tension could be overwhelmed by 
these forces. Finally, the D-brane would become small enough so that the mixing 
amplitude with closed strings is no longer negligible (see appendix |C]). In fact, 
way before this happens, the string would have absorbed 0(N) Z fields so that the 
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planar approximation used in the calculations will be invalidated. The outcome of 
the instability is then difficult to predict but can have important consequences in 
string theory. 

4.1 One Site 

We start with the Hamiltonian for L — 1: 



H = 2\ \a> 2 + a} a + a(a) + a 



(4.1) 



The advantage of this case is that we can diagonalize H explicitly for q = 0. This 
is because a) and a act as shift-left and shift right operators for q = 0. Thus, the 
matrix representation of H looks like the following matrix 



H ~ a 2 + 



/0 a ...\ 

a \ a . . 

a 1 a ' ' 

v '■ '■ '■ '■/ 



(4.2) 



and this is the same hamiltonian for a particle hopping on a semi-infinite lattice, as 
well as the quadratic form that one considers when studying the potential energy of 
a semi-infinite collection of beads attached by springs. 

The complete set of eigenstates consists of the ground state ( |3.11|) with L = 
1, which is considered as a bound state for the hopping particle to scatter of the 
boundary and the continuum, 



l^(^)) = /. [sinfcn + asmk(n + 1)] \n) , with < k < n , 



71=0 



and energy 



E(k) = 2X(l + 2acosk + a 2 ) 
The gap with the ground state turns out to be 

2A(l-a) 2 . 



(4.3) 



(4.4) 



(4.5) 



Moreover, the maximal energy of the continuum is 2A(1 + a) 2 . These continuum 
states follow a delta-function normalization 



(*(*0I*(*)> 



71 



I + 2a cos k + a 2 )5(k' - k) . 



We can then use them to build normalizable wave-packets 



10) 



dkf(k)\^(k)) 



(4.6) 



(4.7) 
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choosing a function f(k) so that (<j)\4>) = 1 and the initial occupation number 
(0|n(O)|0) is finite. Then we can have normalized states with finite energy but 
with growing occupation number. 
For instance, the wave-packet 



v / dhsmk\V{k)) = |(a|0)H-|l}) , 



(4- 



increases its mean occupation number while it evolves, as it is depicted in Figure |l|. 




Figure 1: Mean occupation of wave-packet (4.8). 

This continuum spectrum feature is not only a property of q = 0. Indeed, when 
q < 1 we always find the same behavior. This is because for sufficiently large n 
(occupation number), the hamiltonian asymptotes exponentially quickly to a similar 
form 



H 




(4.9) 



since q n — > 0. 

4.2 Two Sites 

For two sites we can no longer solve for the exact spectrum but we can do (degenerate) 
perturbation theory around the maximal case (a = and q = 0). The unperturbed 
Hamiltonian for two sites is, 



H = 2A a\ai + 2A a' 2 a2 — A(d|a 2 + aic4) • 



(4.10) 



Since this Hamiltonian correspond to an integrable spin chain, we can diagonalize 
it using the usual Bethe Ansatz. The eigenstates are (e.g. see ||25|]), 



\¥°\k)) = A n (k) J2 [2sin(H) - sm{k{l + 1))] \l,n) 



(4.11) 



1=0 



-25- 



where the states |m, n 2 ) have been re-labeled |ni, ni + 712) = |/, n.) to reflect the fact 
that the number operator commutes with Hq. The energy of these eigenstates is 
given by, 

E(k) = 2A(2-cosfc) , (4.12) 

and the complex momentum k is determined by the equation, 

4 sm(kn) - 4 sin(£;(n + 1)) + sin(k(n + 2)) = . (4.13) 

We now want to consider first order perturbation theory in a around these eigen- 
states. Again, the perturbing potential is ( p.20|) and we ignore the constant term in 
( |3.18|) which is 0(a 2 ). In particular we are interested in the continuum energy bands 



that could arise from turning on the perturbation. For that to occur we will need to 
have an infinite degeneracy or an accumulation point in k in the large n limit. 

Let us start by looking at the infinite degenerate roots of ( |4.13j ). It is easy to 



show that for a given n, the degenerate roots have the form k = air/b with a,b G Z 
and where the degeneracy occurs for states \^ n , ) with n' = n + Z6. Now, the 
perturbation ( |3.20| ) has non-zero matrix elements only between unperturbed states 
differing in n by one. Therefore, at first order in perturbation theory we only need 
to consider the degenerate roots: k = 0,tt (we can take k G [0, n] without loss of 
generality) . 

For k = 0, the unperturbed energy is E^ = 2A. To find the first order correction 
in a we need to diagonalize V in the subspace of the unperturbed energy eigenstates 



( 4.11 ). However for our purposes, we only need to look at the large n limit to see 



the presence of the continuum. One can show that in this limit, 

m?V 2Aa (|*£0 + l*£i>) . (4.14) 

Thus the asymptotic form of the first order correction to the eigenstates is going to be 
of the form of a superposition of plane waves \^n ) ~ ^ n e ipn \^n ) with continuum 
momentum p and energy: 

£ (1) =4a\cosp. (4.15) 

If we follow the same procedure with the root k = n we find that at large n the 
matrix elements of V vanish and so there is no continuum band at this energy. Next, 
we consider the case of complex k and we find an exponential accumulation point 
at k — i log 2 with energy E^ = 3A/2. The asymptotic form of the action of the 
perturbation is, 

mi 0) > - \^ (i*Si> + i*£i>) > ( 4 -!6) 

and so we get a continuum band with, 

£ (1) = 3aAcosp. (4.17) 
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4.3 Multiple Sites 

For multiple sites a quantum analysis is no longer feasible. Nevertheless the presence 
of continuum bands can be seen directly from a classical analysis. For any number of 
sites the classical limit is ft — > and thus we can use the coherent state Hamiltonian. 
For completeness, we will work with < q < 1 and then take the limit q —>■ when 
needed. 

Looking at the action (|3.24j ) we see that our system is subject to the constraints, 



V n = ° . 
P* + /«(rO=0, (4.18) 

and then, the phase space is isomorphic to the configuration space fa, fa), which is 
a product of L discs D q of radii j^-. Using the auxiliary constants r = r^+i 



a 



and (f>o = 4>l+i — ti" the classical Hamiltonian can be simply posed as 

L 



H L = 2\J2 n (r-i - n+i cos(0 4 - <p i+1 )) . (4.19) 



i=0 



In particular, let us first consider the one-site classical Hamiltonian 

Hi = 2A(ft 2 + r 2 - 4ftr cos <p) . (4.20) 



Figure 2: Closed and open orbits in D q . 

Curves of constant energy, which is the unique conserved quantity, determines 
the trajectories in the two-dimensional phase space D q . In the present coordinates, 
the Tli = E curves are nothing but circles whose centers are displaced a distance a 
from the center of the disc D q (see figure 0). The radii of these orbits are given by the 

energy through yjE/2\. As long as E < 2A(^ ft) 2 the orbits are closed circles. 

However, for energies 2A(^^ ft) 2 < E < 2A(^z ft) 2 the orbits are open arcs. In 

a semiclassical approach, energies for which the trajectories are open in the phase 
space give rise to a the continuum in the quantized system. The proper measure of 
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the continuum is the phase space area between energies E and E + 5E. This area 
has to become infinite once the trajectories reach the boundary (otherwise one would 
find that the phase space has finite area, and then the total number of states would 
have to be finite). 

The minimal and maximal energies of the continuum in the q — > limit are seen 
to agree with the results of the quantum perturbation theory. 

We now finally come to the case of arbitrary number of sites. In this case we 
will only calculate the energy gap between the ground state and the first continuum 
band. To obtain the energy for which the continuum begins, one has to look for the 
minimal energy for which the energy constant hypersurfaces intersect the boundary 
of the phase space D^. To this end one needs to compute and compare the absolute 
minimum of the resulting L functions when one of the r« is set to -^— . These L 
functions are bounded from below by quadratic polynomials on the radii r^ 

L 

K = HlL = i > Mi = 2Aj> (r, - r w )| i , (4.21) 

J 1 — q ( * J 1— q 

j=0 

and coincide with them when all angles fa are taken to be n. Then it suffices to 
find out the minima of the Ai 3 L . The stationary points should satisfy in all cases the 
recursive equation 

in - r i+l - n-i = , (4.22) 

subject to the boundary conditions 

r Q = a, rj = - , r L+1 = a. (4.23) 

1 — q 

The recursive equation (|4.22|) is fulfilled by a linear dependence of r^ on the site 
labeling i, with the constants adjusted to also fulfill the boundary conditions ( |4.23| ). 
The critical point of each A4 J L is 

l-a(l-g)\ 

i + a if % < j 



j(i -q) 

a(l-q)-l \. L+l-ja(l-q) 



(L+l-j)(l-q)J (L+l-j)(l-q) 



if i > j . (4.24) 



It is easy to see that each critical point is a minimum of each A4 J L , since all the 
eigenvalues of the corresponding Hessians are strictly positive. The energy evaluated 
in these points is 

H* I = KL + l)(l-a{l-g)f , , 

As a function of j, that takes values 1 < j < L, the minimal energy is obtained for 
the j corresponding to the central site. So, the minimal energy for which constant 
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energy hypersurfaces intersect the phase space boundary is 



E, 



cont 



4A(l~q(l-g)) 2 

(L+i)(i- q y 

4A(L+l)(l-q(l-< ? )) 2 
L{L+2)(l-qY 



if L is odd 



if L is even 



(4.26) 



Let us now turn back to case q — > 0. For L = 1 we have again E cont = 2A(1 — a) 2 , 
which is the minimum of the unique band ( |4.4j ) of the quantum system. For L = 2, 



we have E, 



cont 



3A/2(1 



a } 



which is the minimum of the lowest band ( [4.171 ) 



computed perturbatively for small a. We are eventually interested in the limit of 
a large number of sites L. According to the classical analysis, the minimum of the 



lowest band is, in that case, 



E r 



4A(l-a) 



(4.27) 



As we argued, from this energy the quantized system will have a continuum spectrum. 
Then, this is the amount of energy needed for having configurations with a mean 
occupation number growing in time. Note that this energy is accessible to long 
semiclassical strings. 

Let us end this section, presenting how the transition from oscillatory to non- 
oscillatory solutions near the value E cont is observed in numerical solutions. In figures 
H, |] and H we show different classical solutions of the system with L = 30 and their 
corresponding occupation numbers. The 3 cases were solved with similar initial data. 
Only the amplitude of the initial profiles (r i; fa) were slightly varied to change the 
energies. Moreover we chose the energies to be near the critical value E cont . For 
energies below E cont , no ti can tend to 1, so the mean occupation number is kept 
bounded. But as long as the energy is bigger than -E con t it is possible to find solutions 
with some r^ approaching 1 and a mean occupation number growing monotonically 
in time. 



<n> 




Figure 3: L = 30 solution with energy E = 0.94.E, 



cont' 



4.4 A toy model for the D-brane instability 

We have found, surprisingly, that the spin chain model associated to a string attached 
to a non-maximal giant exhibits a continuous spectrum. Naively, this is counter to 
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<n> 



Figure 4: L = 30 solution with energy E = 0.99-E, 



cont • 




<n> 




Figure 5: L = 30 solution with energy E = 1.05-E, 



cont- 



our intuition from AdS, as we have a finite energy state, and in the dual field theory, 
there are only finitely many states with energy less than E. This means that the 
true spectrum of the string plus D-brane system ultimately must be discrete. 

Considering how we got our spin chain model, we have to notice that we had 
various assumptions built in our computations. First, the giant has energy of order 
N, while the string should have an energy of order y/N, roughly, this corresponds to 
occupation numbers lower than y/N. This is required for the planar approximation 
to Feynman diagrams to make sense. If the total angular momentum of the spin 
chain becomes large with respect to yN, the planar approximation breaks down. 
This means that the Hamiltonian in terms of the boson spin chain we have described 
is only valid in some regime (where the occupation numbers are not too large), and 
when we get out of that regime, the fact that N is finite is important. As we take N 
large, the bosonic spin chain model realm of applicability grows. So the continuum 
spectrum and bands we have been discussing are a property of the strict large N 
limit. At finite N, the continuum spectrum is resolved into a discrete spectrum 
with an extremely fine spacing of the eigenvalues. This would be very hard for an 
observer to determine. This continuum spectrum, which is associated with large 
quantum numbers for the occupation numbers of the chain has to be interpreted in 
terms of classical physics (this is the same reasoning that went into the work ||54j| ). 
At large quantum numbers, the problem becomes classical, and we should be able to 
come up with a classical argument that shows that the strings should grow in size, 
so long as string interactions and the brane back reaction can be ignored (these are 
the 1/N effects that the planar approximation does away with). 

Now, we want to come up with a model of the system that clearly shows the 
instability, without the complications of the AdS^ x S 5 geometry. The idea is simple: 
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we have a charged D-brane in the presence of a magnetic RR background. Because 
of the background, a moving D-brane is accelerated with respect to geodesic free-fall 
by a Lorentz force. Since fundamental strings are not charged with respect to the 
RR-background, they should prefer to follow geodesies. 

This motivates the following toy model for the giant graviton setup: the ends of 
the strings are attached to the D-brane, so they are accelerating. By the equivalence 
principle, we can think of the same system as a string with fixed ends suspended in 
a gravitational field. A close analogy of the setup is the act of suspending cables 
between telephone poles. The telephone poles are replaced by the D-brane itself , and 
it is supported in place by the RR-background. The cable between the telephone 
posts is the string. Unlike conventional cables, the string tension is constant (a 
fundamental constant in perturbative string theory). This means that in principle 
the string can weigh more that the tension of the string can support. This requires a 
long string (long enough so that the weight of string between the two posts is bigger 
than the string tension), whose length depends on the strength of the gravitational 
field. 

Since the string can not break (we are assuming that we are in the strict planar 
approximation), the string will grow by falling, if it is long enough to begin with 
(this can be accomplished by separating the putative telephone posts sufficiently). 
In fact, the string will keep on falling forever if string interactions are turned off. 
Otherwise, the string will get so long that eventually there is a finite probability per 
unit time that the one long string will break into a closed loop plus a shorter string 
suspended between the post (this can be understood as the emission of gravitational 
waves, or other fields, by an accelerated object). 

Our toy model clearly exhibits the instability of the string system that we found 
with our spin chain model. It also provides for the presence of the gap: if the 
telephone posts are close together, the minimal energy string will be stable. Indeed, 
many small perturbations of the configuration will be allowed. To make the string 
fall, we might need to stretch it a lot, but once we have enough weight in the system, 
the string can not sustain the total weight of the string, and the instability we discuss 
will have set in. Because we have a regime where the string is stable, this regime 
of lower energy will have a discrete spectrum, and when we reach the instability, we 
will get a continuous spectrum. This gap is determined by the distance between the 
telephone posts (this can be related to the total angular momentum of the string 
ground state), and by the effective gravitational field. This effective gravitational 
field is determined by how much is the giant graviton accelerating: the bigger giants 
move slowly, so their effective gravity is small, while small giants move fast and 
accelerate a lot more. 

This means that the gap in the spectrum should become smaller as the D-brane 
becomes smaller (in our spin chain notation, a becomes bigger). This is exactly what 
we see from eq. [4.5|, as well as from understanding the details of figure 0. 
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5. Evidence for Integrability 

In this section we present numerical evidence for the integrability of the classical 
Hamiltonian ( f4.19|) for q G [0,1]. For L sites, this would correspond to have L 
constants of motion 5 . Then, trajectories in the 2L-dimensional phase space take 
place in a L-dimensional torus. However, with the exception of separable systems, 
there is no systematic procedure for finding constants of motion. Even in the rather 
simple case of two sites, we do not know another constant of motion apart from the 
energy whenever a^O. Hypersurfaces of constant energy are 3-spheres in D^, and 
if a second constant of motion existed, the motion would be in a 2-torus included 
in the 3-sphere. Then, the intersection of this 2-torus with a any hypersurface (of 
dimension 3) would be in general a closed curve in the phase space. We can verify 
if this is the case by studying numerical solutions of the system. For example, in 
Figure |6] we plot the values of coordinates (r 2 ,0 1 ,0 1 ) of a given solution for many 
different times when the variable r\ takes the value 1/2. The fact that points lie in 
a closed curve is a strong indication that the motion is taking place in a 2-torus and 
the system is integrable. One can try this for many different hypersurfaces and one 
always gets closed curves. Similar evidence is found for cases with higher number of 
sites. 




Figure 6: Intersections of the a solution with r\ = 1/2. 



Note that for q = 1 the deformed operators (E.l|) become the usual harmonic 
oscillators and the system becomes a lattice of ordinary bosons. In this case quantum 
and classical integrability follow trivially. For q ^ 1 the classical Hamiltonian ( [4.19| ) 
can be seen as a deformation of the Hamiltonians studied in |55, [56], 57fl . We now 
conjecture that our family of Hamiltonians parameterized by q, a is indeed a two 
parameter family of (quantum) integrable Hamiltonians. 

Because the spectrum of the quantum hamiltonian has a continuum, it is very 
unlikely that the system can be solved via a simple Bethe Ansatz for all a, even 
though for a = the system can be solved via a Bethe Ansatz. The distinction 



5 Moreover the constants of motion should be compatible with each other, which means that the 
Poisson brackets with each other should vanish. 
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between a = and a 7^ resides in the fact that the constant energy surfaces 
are compact in the first case, while they are non-compact in the second case (they 
intersect the boundary). This non-compactness of the energy surfaces for the spin 
model upsets the usual Bethe Ansatz intuition for finite systems. In finite systems 
the Bethe ansatz hinges in having some discrete conserved number of quasi-particles, 
and this produces a system of algebraic equations that describe scattering of these 
quasi-particles with each other and the boundary. The solutions of these Bethe 
equations give a finite number of roots that one identifies with the discrete spectrum 
of the finite system. 

To have a continuous spectrum one needs one of two conditions: either the 
spin chain is infinite, or the conserved "number of quasiparticles" is described by 
a continuous parameter. The first case does not describe our system well, as we 
have found that the ground state has a well defined finite length (average number 
of spins). The second case would not produce algebraic equations that solve for the 
spectrum and we would be hard pressed to call such a system a Bethe Ansatz. 

However we seem to have an integrable system with a continuous spectrum. We 
should contrast this observation with the observation of [27| that the Bethe Ansatz 



seems to break down at two loops for the maximal giant graviton, or that the BMN 



limit breaks down [28]. There are also other hints that the asymptotic Bethe Ansatz 



conjecture is incomplete [|3~T |. We interpret all of these facts not as a breakdown of 



integr ability, but instead as a breakdown of applicability of a Bethe Ansatz. 



6. Discussion 

In this paper we have shown an example that dynamical finite D-branes can be 
treated consistently in the AdS/CFT setup. In particular we showed that we were 
able to characterize the D-brane in the most conservative form: as a geometric locus 
where open strings can end. We saw that this property could be derived exactly 
from the dual field theory. We were able to do this not just for a maximal giant 
graviton, but also for smaller giant gravitons that move in the AdS spacetime. The 
calculations required in field theory were rather subtle and we found a lot of surprises 
in trying to explain the dynamics of this system in detail. 

The first surprise we found was that the trajectory of the spherical giants we stud- 
ied contained families of null geodesies, and that it was possible to take a standard 
Penrose limit along these trajectories. We found that this limit gave us a different 
coordinate system than the usual one, because it was adapted to the D-brane being 
static. 

We then found that is was possible in the limit of large quantum numbers to 
get the classical ratio of different components of angular momentum carried by the 
ground state of the string. Part of it was along the direction of angular momentum of 
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the giant, and the string and giant in principle exchange angular momentum between 
them. 

We then looked for a description of these semiclassical states of the AdS geometry 
directly in the dual field theory in the perturbative weak coupling approximation. 
We found a spin chain model with variable numbers of sites at the boundary. We 
were also able to describe this system in terms of a boson chain model, where the 
individual bosons obeyed Cuntz algebra oscillator relations. The Hamiltonian we 
found was quadratic in generalized raising and lowering operators. It results as a 
limit q — » of a g-deformed chain of harmonic oscillators with nearest neighbor 
interactions. The boson system did not preserve the total boson occupation number, 
so the boundary conditions we found are of non-diagonal type. 

A big surprise was that when we tried to diagonalize the boson chain model we 
found a continuous spectrum of the effective hamiltonian. This was very unexpected 
when we started the project, but we found an analogy that made it intuitively obvious 
and necessary for the AdS/CFT correspondence to work. The analogy consists on 
thinking of the system as an open string whose ends are accelerated because they 
are attached to a D-brane, in the presence of a RR background. The background 
accelerates the D-brane, and the D-brane drags the string itself. The continuous 
spectrum was argued to be related to an instability that appears only if the string is 
too long and can not support its own weight. This instability might be interesting 
for the study of cosmic strings, because it can serve as a mechanism to seed long 
strings in an expanding universe. 

We also found numerical evidence that the open string attached to the giant 
graviton was in general to one loop order governed by an integrable system. Our 
evidence is numerical and classical. We also found that the continuous spectrum 
seems to prevent the integrability of the system to be described by a Bethe Ansatz. 

Nevertheless, despite this progress, there are many loose ends and open problems 
left to study. Here we would like to list some of them to encourage future research 
along these lines. 

1. The Plane Wave Spectrum and Integrability 

The issue of the integrability of the Hamiltonian ( |3.9| ) is important to prove 
the presence of the BMN spectrum. However, the fact that the spectrum 
has continuum energy bands suggests that integrability might be realized in 
a way quite different from the usual Heisenberg spin chain. In fact, since we 
have absorption and emission from the boundaries, one wonders if the usual 
formalism using the boundary reflection matrices could be generalized to take 
into account this effect [^|. Whatever the answer, this problem can have 



implications beyond the study of the AdS/CFT correspondence and should be 
very interesting to resolve. 
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2. The D-brane Instability 

As we mentioned before, the end state of the unstable giant graviton is difficult 
to predict from our current results. The reason is that all the approximations 
used in our field theory calculations are invalidated when the word's length 
describing the string is large with respect to yN . In any case, the end state of 
this instability can be important for the study of certain classical spacetimes. 
In particular we know that the 1/2 BPS supergravity solutions found in |H| can 
be interpreted as coherent states of these giant gravitons (and the AdS giants). 
One can excite massive open string modes on these backgrounds and ask what 
effect the instability studied here can have on these backgrounds (one might 
imagine that these might serve as a mechanism to transport charge between 
droplets for example). This can be important for the stability of black holes 
since one can imagine that these will be constructed by exciting many open 
string modes on the D-branes. Finally, one can ask if this instability is also 
present in the case of AdS giant gravitons. 

3. Beyond the SU(2) Sector 

In this paper we only considered open strings with two angular momenta on the 
sphere. In other words, an SU{2) subsector of the corresponding gauge theory 
operators. However, it is interesting to go beyond this sector. At this moment 
we do not know how to write down a "nice" Hamiltonian describing the full 
£0(6) excitations on the sphere. However, one can consider an SU(3) subsector 
by using the three holomorphic fields: X, Y and Z. The generalization of our 
Hamiltonian is straightforward. The most general word can be labeled as, 

(YZ^XZ^Y ■ ■ ■ Z nL X)j = \s u m, s 2} n 2 , . . . , n L , s L+1 ) , (6.1) 

where Si =| or J, = X orY, is an SU(2) spin label and in is the usual bosonic 
occupation number. 

The interactions between the sites are much the same as before except when any 
rii = 0. In this case we have the additional permutation interaction between 
the two spins s« and Si+i (see Eq. flC.14|) ). One can then write the Hamiltonian 

as . 



H = H hoson + \J2( 1 ~ Am)(l - ofoj) , (6.2) 

i=i 

where H\> oson is the Hamiltonian for the bosonic lattice ( |3.9| ), and Pij,+\ = 
2(1 + 45; • Si+i) is the permutation operator acting on the spin sites. 
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Constructing the action for the coherent states goes as usual (see appendix [D]) 
We get, 



S = -L I dt J da 
A 



r 2 • 1 A 



+ - cos 6 (p + —(r u + r z (j)' 2 ) 



\ -r 2T 2 r L 2 

(6.3) 



+ ^{l-r 2 ){6' 2 + sm 2 6f' 2 

where the angles 9, if label the SU(2) coherent states. On the other hand, 
one can study the dual string theory and write the Polyakov action for open 
strings as we did in section ||. Using the natural coordinates adapted to the 
brane, and making the expansion in large total angular momentum in the X 
and Y directions we get the action (|6~3|) with the exception that the last term 
has a different r dependence: ~ J dtda(l — r 2 ) 2 {6' 2 + sm 2 9(p' 2 ). It would be 
interesting to understand the nature of this discrepancy. 

4. Multiple D-branes 

Recently it was shown how to write down the operators corresponding to multi- 
ple giant gravitons with strings attached to them |23| . In general the combina- 
torics are hopelessly complicated. However one can study the case of two giant 
gravitons with two strings stretching between them (we need two strings to 
satisfy the Gauss constraint [E3|). Intuitively one would expect Hamiltonians 



like ( |3.9| ) but with more general coefficients for the boundary terms (different 
values at each boundary). In this case there can be a "current" of Z fields from 
one giant to the other. It would be interesting to see what happens in this case. 
We expect that the giants attract each other. We would also expect the same 
instability we found for these models. Moreover, one should be able to measure 
the distance between the two giants using the sigma model representation of 
the coherent states. More generally, one should also be able to study D-branes 
at angles, and measure the angles between them by finding the spectrum of 
open strings stretching between them. 

For many giants the easiest description seems to be in terms of the matrix 
model of j|I|. However it is still not clear how to introduce the open string 
excitations in that language 6 . It would be interesting to re-derive the results 
presented in this paper using the matrix model formalism. This can also give 
a better understanding of the backreaction of the giant gravitons with massive 
open string modes. 

5. AdS Giants 

One can try to extend the results of this paper to the study of the giant 
gravitons that expand in the AdS space. The operators are similar to the 



D See pq for the emergence of closed string BMN excitations in the matrix model 
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ones considered here but with the e symbol replaced with totally symmetric 
tensor contractions. The combinatorics can be done in a similar way but the 
field theory interactions are more complicated when we consider an open string 
spinning along the AdS directions. This is because, instead of scalar "letters" 
we need to use covariant derivatives. However, one expects a simple description 
in terms of spin chains (or a bosonic lattice) as for closed strings |5l], |TD|. 

6. Giants of critical angular momentum 

In this paper we have considered giant gravitons whose angular momentum 
grows proportionally to N. It would be interesting to study gravitons with 
angular momentum p ~ \/N. For p 2 /N fixed and p 2 /N <C 1 the gravitons 
should be described as point-like supergravity modes (as it is done in the BMN 
limit). However, for p 2 /N fixed and p 2 /N ^> 1 the appropriate description is 
in terms of expanded D3-branes. From (|2.6|) one immediately realizes that the 
angle 9$ goes now to zero in the large N limit and each element of the giant 
is traveling in an almost null geodesic. Moreover, from (|2.7|) it is evident that 
the radius of the giant remains finite in the limit. Then, a short open string 
attached to this giant is effectively described by an open string attached to a 



spherical D3-brane of finite radius in a pp-wave background |6T |. The string 
does not need to be spinning along the giant. Then, bosonic lattices with a 
finite and small number of sites should be used to describe the dilation operator 
in the dual field theory. 
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A. Penrose limit 

The Penrose limit is accomplished defining new coordinates by means of a linear 
transformation and taking an appropriate R — > oo limit. The coordinate playing the 
role of the curve parameter, must appear with the same coefficient in t, ip and 77, 



if we want to capture the geometry near the null trajectory ( |2.9| ). The rest of the 



coefficients can be fixed by demanding the R —>■ 00 limit to be well-defined. Then, 
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we consider the transformation 



t = u+- 2 , 


P= R 


V X 

1p = U~ D^ +taIl6l 0-D , 


e = e 0+ y- 


V X 


z 


R sin On 



(A.1) 



After this scaling and the limit R — > oo, the metric becomes, 



ds 2 = —Adudv + Aydudx — (r 2 + z 2 )du 2 

+dx 2 + dy 2 + dz 2 + z 2 d£ 2 + dr 2 + r 2 dtt 2 3 . (A.2) 



If we define cartesian coordinates 



Z\ = z sin £ , z 4 = r sin yj cos r\ 

z 2 = z cos £ , Z5 = r cos if' sin £' , 

z 3 = r simp' sin r/' , z 6 = r cos <p' cos £' , (A. 3) 



the metric is written as 



ds 2 = —Adudv + Aydudx — 2_, z 1du 2 + dx 2 + dy 2 + N^ dz 2 , (A. 4) 



a=l a=l 



and the RR 5-form field strength is 

F(5) = 2dw A (dzi A c?z 2 A dz 3 A c?z 4 + cbs A dz$ A c/z 7 A dz$) . (A. 5) 

This pp-wave configuration can be put in the standard form using the following 



coordinate transformation 37 



x + = u , x\ — x cos u + y sin u , 

1 
x_ = v — —xy , X2 = —x sin w + y cos it , 

x a +2 = z a , for o = l,..., 6. (A. 6) 

which leads the metric to 

8 8 

ds 2 = —Adx + dx- — y^ x 2 dx\ + 2_, dx 2 . (A. 7) 

Notice that this transformation involves rotations at constant angular velocity in 
the x, y plane. In this sense, the D-brane we will be considering are rotating in the 
standard plane wave limit. 
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B. Brane in the Penrose limit 



To see how the D3-brane is specified in terms of the new coordinates ( |A.1| ) we should 
apply the coordinate transformation to the boundary conditions and then take the 
R —>■ oo limit. In the original coordinates, the giant expands in ((p,rj,£,), and then 
for an open string, they should satisfy Neumann boundary conditions 



ty o-<,| cr= o, 7r 











On the other hand, the giant is situated in p = and 9 
coordinates have Dirichlet boundary conditions 






0. 
0. 



(B.l) 
(B.2) 
(B.3) 

8n. Therefore, these 



(B.4) 
(B.5) 



The remaining transverse direction to the D-brane is ip and the D-brane is moving 
along this angle with ip — t. The corresponding Dirichlet boundary condition is 

*VW = 6t l=o,« ■ ( B - 6 ) 

As a consequence, the remaining Neumann boundary condition also mixes t and ip 



d a t\ a= ^ = cos 2 6q d^\ u= ^ 



(B.7) 



Using (|A.1|) the set of boundary conditions ( |B.1|) - (|B.7| ) is translated, before the limit 
R — > oo, into 



^o-Z\ cr=0 







ftrUL n _ = , 

u l(7=U,7r 

, cot O 

a ff ul _ nu = 

u ltJ=U,7r 



i? 



<9 CT x 



cr=0,7r ' 



lcr=0,7r 



<^|<r=0 )7 r 



2 cot O 



5i> 



<T=0,7T 



(B.8) 

(B.9) 

(B.10) 

(B.H) 
(B.12) 
(B.13) 

(B.14) 



So, in the R — > oo limit, we have Neumann boundary conditions for w, ?;, Zi, 2^ and 
Dirichlet boundary conditions for x, y, z 3 , z±, z 5 , z 6 . 
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C. Combinatorics and field theory calculations 

We begin by listing some properties of the totally anti-symmetric tensor that are 
useful for field theory calculations. This tensor is defined as 

{1 if (ii • • • i p ) is an even permutation of (j\ • • ■ j p ) 
— 1 if {i\ ■ ■ -ip) is an odd permutation of (ji • • ■ j p ) , (C.l) 

otherwise 

where p is any integer p < N, and i\, • • • ,i p and ji, • • • ,j p are integers from 1 to N. 
The simplest examples are 

i- = 5 i , 

4 = *W - &A ■ (c-2) 

Some of most useful properties of the e tensor are 



31- jp 


/ -A / Jx 3l---3x-i3x+\---] P ' 
x— 1 


(C.3) 


h---ikjk+i---j p 


_ (iV-p + fc)! lfc+1 ... ip 


(C.4) 


Ji—ih Ji-Jp 
ji—jhh-lp 


_ f £ ] e i l---i k jk + l---3p 
h ... i p ' 


(C.5) 



Another useful identity that can be derived from the ones above is 

h-jp-i-yi n-i P -i72-7fc _ (g ~ l)-(-^ — fc + !)• / . 7l l2 - lk _ {P ~ 1) 7 i- 7fc 
ii—ip_iw ji—jp-iiia—Pk fjy _ „ _ ^ i 2)1 I i" 1 M2--M*: (jy _ £ + 1") Mi-'-Mfc 



(C.6) 



In the rest of this appendix we present some of the field theory calculations 
leading to the Hamiltonian ( p.9|) . Most of the calculations are very similar to those 
presented in |25|] in the context of maximal giant gravitons. 

We begin by pointing out that an operator with a Z field at the end or begin- 
ning of the word does not represent a linearly independent state. In fact, from the 
identities above one can show that 

<X Z l ■ ■ ■ Z t-K ZW t = \<-Z Z n ■■■ Z t ^(W) - \<-tl Z h ■ ■ ■ Z l W tl ■ 

(C.7) 

This shows that the operator with a Z at the border of the word factorizes into a 
brane with a closed string and a bigger brane. 

Now we turn our attention to the correlation functions. These will be computed 
to one loop in the 't Hooft coupling and using the planar approximation of the large 
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N limit. We will be using the bosonic part of the SYM action given by 

S = ^- [ d 4 xTi {^F^F^ + D^XD^X + D^YD^Y + D^ZD^Z + V D + V F J , 

(C.8) 
where, 

Vb = ^Tr (\[X,X) + [Y,Y] + [Z,Zf) , (C.9) 

^ = 2Tr(|[Y,Y]| 2 + |[X,Z]| 2 + |[Y,Z]| 2 ) . (CIO) 

Let us begin with the free field correlation functions, since they are going to give 
us the normalization of the operators. The propagator for the scalars is of the form 

<#(s)#(0)> = ^lA-M > ( c - n ) 

J Z7T |X| Z 

where is any of the complex scalars X, Y, Z. Thus we see that the correlation 
function of any word with a total classical dimension Ao will be multiplied by the 
following overall numerical factor: (J^ T^n)^ - In what follows we will drop this factor 
from the calculations for simplicity. 

Let us now compute the correlation function for the operator 

where the word W has classical dimension L. We will assume that the operator is 
very long such that p > yN and L < yN with p ^> L. In the large iV limit, the 
leading contribution will come from contracting all the Z fields and then the words 
planarly. Moreover, in the large iV limit operators with different values of p and 
word lengths will be orthogonal. Thus the free theory two point function gives 



\ W W W W/free \P L )- t ii-i p -iHi t ji-jp-iM v K 7i VV 72 /free 



(p-l)\ 2 (N-2) 

(N-p-l)\ 

p\(p - 1)\(N - 1)\N L 
(N-p)\ 



(P-1) 

N -p 



(Tr(W) Tr(W)) {ree + K -^-(Tr(WW)) iree 



(C.13) 



In going to the second line of (|C.13|) we used the identity (|C.6| ) and for the third line 



we contracted the words using the planar approximation. 

If the word W is made of a single type of field we will get an additional multi- 
plicative factor of L(l + 0(L/N)) to ( |C.13| ) from the cyclic property of the trace. 
This is one of the reasons we consider only L <^ N and in particular L < yN . 

We now turn on the Yang-Mills coupling and consider the interactions between 
the scalars. Let us start by calculating the interactions in the bulk of the words. This 
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terms are obtained when all the Z fields of the giants are contracted we free field 
propagators and the vertex is entirely contracted with letters of the words. In the 
planar approximation these will be the familiar nearest-neighbor interactions that 
were considered by Minahan and Zarembo ||. Thus if we regard each "letter" in W 
as a SO (6) vector, the nearest-neighbor interactions have the familiar form, 



H 



1 



i+i 



\[K U+1 + 2(/ M+1 - P M+1 )] , 



(C.14) 



where K and P are the trace and permutation operators respectively and A = 
g s N/2n. For the SU(2) and 577(3) sectors considered above, the trace will be zero 
and it is easy to see that the remaining interactions account for the first two terms 
in the Hamiltonian (|3.9|) . 

For completeness, let us show how to calculate the bulk part of the Hamiltonian 
by considering the example of the interactions giving rise to the identity in ( |C.14 ). 
To this end we study the interaction between a X and Y fields of the operator 



rnv — Ji—3p 7*i 



Zf-_ 1 i {W 1 XYW 2 



jp 



(C.15) 



It is well known that for chiral operators the interactions coming from Vd cancel with 
the gauge boson exchange and the scalar self energies. Therefore, the only relevant 
interactions between X and Y will come from the F term. It is not difficult to see 
that the only planar interaction will be, 



/ d 4 x Tr (XYYX) . 

Kg s J 

This interaction is illustrated in figure [7|. 



(C.16) 



Wi x 



W: 




W, X Y W 2 



Figure 7: One-loop interaction between two nearest-neighbor letters of a word generated 
by a vertex of the form Tr (XYYX). 



The rest of the free theory contractions go as before. Therefore, all we get from 
the interaction is a multiplicative correction to (|C.13|) of the form 



g s \~ 2 (9. ^ 



2n 



, N 
2ixJ \ %g s 



x\ 4 / d 4 y- 



y\ 4 \x — y\ A 



-2Alog(|x|A). (C.17) 
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Thus, after dividing by the norm ( |C.13|) and using the definition of the anomalous 
dimension matrix (|3.5| ), we see that this interaction will give a numerical contribution 
of A to the identity interaction in agreement with ( |C.14| ). 

Now let us see what happens when we have a non-chiral operator. Let us consider 
the word W\XYW2- In this case the D-term interaction will no longer cancel with 
the gauge boson exchange and the scalar self-energies. However, we can exploit the 
fact that these last two interactions are flavor blind to write their contribution in 
terms of the D-term contribution to the chiral operator. This is illustrated in figure 



= 




x 




gluon + 
self energy 



X 
X 



Y 
Y 



gluon + 
self energy 



X 



Y 



Figure 8: The first line of the figure illustrates the cancelation of the D-term contribution 
with the gluon exchange and the scalars self energies for a chiral operator. Since the gluon 
and self energy interactions are flavor blind they will give the same numerical value if we 
compute them using the non-chiral operator as in the second line. 

It is not difficult to see that there is no planar contribution to the identity coming 
from the F-term. Thus, the total contribution to this interaction is from the diagrams 
shown in figure || 



x 




x 



x 




Figure 9: Total contribution to the one-loop interaction between X and Y of the word 
W\XYW2- The second term is the gluon and self energy interactions calculated in terms 
of Vd using the result illustrated in the previous figure. 
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The first diagram comes from the vertex 

— Tr(XYYX) , (C.18) 

2rrg s 

and the second from 



Ti(XYYX) . (C.19) 

Note that the second vertex is to be evaluated between the chiral operator. At the 
end, we see that we will get a numerical contribution of A/2 from each interaction, 



and this will sum up to the contribution shown in (|C.14|) . The rest of the interactions 



in the Hamiltonian (|C.14j ) can be obtained in a similar way. 



Let us now turn our attention to the interactions between the word and the Z 
fields of the giant graviton. The first contribution comes from the scalars at the ends 
of the word. Thus we can study the operator 

Ol = <Z Z n--- Z V,^ W Z^ ( C - 2 °) 

where is any of the complex scalars except Z and the word W has length L — 1. 
Doing the free contractions we get 

x(^(x)^(O)[0(O)0(x)]^). 

(C.21) 

We now note that there are only two relevant planar interactions. These will 
have the general form: 

° Tr(ZZ00) , — — Ti(ZZ<f)0) . (C.22) 



ng s Trg s 

If the field is X, Y or Z, the constants a, (3 will need to include the contribu- 
tion from the gauge boson exchange and the scalar self energies. In any case, the 
interactions ( |C.22| ) will give a contribution to ( |C.21| ) of 



WW ~ -2AiV ^ !(p ( ^, 1)! [^ + p (l - | )] log(MA) , 

(C.23) 



up to multiplicative corrections of order 1/p and 1/N. Thus after dividing by 
the norm ( |C.13| ) we get the following contribution to the anomalous dimension: 
A [ap/N + j3 (1 — p/N)]. Calculating the constants a and /3 is just a matter of count- 
ing how many terms like ( |C22| ) we can find in Vd and Vp. At the end we get that 



for the fields X, Y, X or Y, a = and (3 = 1. Thus the anomalous dimension gives 
A(l — p/N). Translating to the boson lattice language and taking into account the 
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contribution to the field at the other end of the word, this gives the first term in 
the second line of the Hamiltonian ( |3.9|) . On the other hand, if = Z we have 
that a = — 1 and (3 = 3/2. Thus the contribution to the anomalous dimension is 
A(l + (l-p/]V)/2). 

Finally, there is an interaction that is not present in the case of the maximal 
giant graviton. This is the exchange of a Z field between the word and the giant 
graviton. In the planar approximation this can only happen when the Z is the second 
or next to last letter of the word. To calculate this amplitude we can consider the 
correlation between the operators 

QP+ i = ^;;iQzi ■ ■ ■ Zf p (XW)l + + \ » ° P = <-'Z Z n ■ ■ ■ Z t-\ ( XZ W)Z , (C.24) 

and assume that the word W has classical dimension L. The calculation here is very 
similar to the previous example. The free contractions give 

<(T +1 (*)(T(0)> ~ N^pl&fc^e^ZiZISiXZX)") . 

(C.25) 

Since these are chiral operators, the only interactions will come from Vp. It is easy 
to see that the only planar interaction will be 

+ — Tt(ZXZX) . (C.26) 

Kg s 

Doing the contractions we get 

where we have used the identities ( |C.4| ) and ( |U.5| ) . 
The norms are 



(V P V P ) iree - K — J — , (C28) 



p\ 2 (N-l)\N L +2 
p(N -p)\ 



Therefore, for large p the contribution to the anomalous dimension is: Ayl — p/N. 
We recognize this as sources/sinks in the Hamiltonian ( p.9|) . 

To conclude we show that closed string emission/absorption is suppressed in the 
large N limit if we take p ~ jN. For definitiveness consider the operators 

O? = €t Z jl ■■■ Z l MY L ) , Ol = 4i:t p Z}l ■ ■ ■ Z%-_\ (YZY L -% . (C.30) 
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Then, the interacting correlation function gives (up to signs or numerical factors) 

Lri 2 (N- IVN L 
(0{Ol) ~ A P ' [ {N _J ] log(|x|A) , (C.31) 

The norms are, 



(Old) 



free 



Lp\ 2 N\N L 
(N-p)\ 



, (0 2 O 2 ) 



P 



2 N\N L 



free 



p(iV-p)! 



Therefore, the contribution to the anomalous dimension is of the order 



A ~ X\ — 

N \ N 



This is clearly suppressed in the limit we are considering of L < yN. 



(C.32) 



(C.33) 



D. Calculations in the SU(S) sector 

Here we give some more details about the calculation of the effective sigma model 
action for the SU(3) sector. We begin with the gauge theory side. The coherent 
states for states of the form fl6.1| ) take the product form 

L 

\CS)=(g)\n t ,z l )®\n L+1 ) , (D.l) 

4=1 

where \zi) are the bosonic coherent states ( |E.2| ) and \Hi) are the SU(2) coherent 

(D.2) 



states defined by [53|, |32 



\n) = e iS ^e iS y e \ |) , 
where n = (sin 6 cos (p, sin 6 sin <p, cos 6) is a unit vector in S 2 . They obey 

1 



(n\S\n) = -n , 
\ i i / 2 . 



(D.3) 



(n\n') 



ros 1(6 - 9') cos -{if -<p')-i cos he + 6') sin -(0 - 5') 



2 V ' 2 

The action for the coherent states is 



e *§(v-^) 



5= lim /dt(i(CS|^|CS)-(CS|H|CS)) , 



(D.4) 



(D.5) 



and we get (|6.3| ) 
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On the other hand, one can study the dual string theory and write the Polyakov 
action for open strings as we did in section 2. In the case of three spins, the natural 
coordinates on the sphere |X| 2 + |y| 2 + |Z| 2 = l are 



X = ±Vl-r 2 cos6e iipi , Y = ±Vl - r 2 sm6 e^ 2 , Z = re*-® , (D.6) 



where the giant graviton is located at r = a/1 — p/N and cp — ip\ + cp 2 — t. Following 
[34}j we can make another change of coordinates: 



1 1 

a = n(Vi + ^2) , = -(<Pi -<p*) , U x = cos 9e lf3 



Uo = sin 1 



-ip 



(D.7) 



In these coordinates, the metric on 1 x S 5 takes the form 



ds 2 



(l-r 2 )dt 2 +- -dr 2 +r 2 d(j) 2 +2r 2 d(j)dt+(l-r 



1 — r 



(Da) 2 + -driidrii 



(D.8) 



where coordinates on a S" 2 are given by rii = £*crj£, where £ = (L^e 10 , U 2 e ta ) and ctj 
are the Pauli matrices. The covariant derivative is defined as Da = da + C, where 
C = -iUidUi. 

We now want to expand the Polyakov action in the limit of large angular mo- 
mentum in a, and choose a gauge where this angular momentum is uniformly dis- 
tributed along the string. Note that this angular momentum is dual to the total 



number of Xs and Ys in the word (|6.1| ). We choose the coordinates in the Polyakov 
action such that a is coupled through its covariant derivatives ~ g D a aDba, where 
D a a = d a a - iUid a Ui. 



With this choice, the lagrangian density in momentum space (2.37) becomes 



C = p t i + p r r + p^cj) + p a D t a + pihi 



(D.9) 



and we have the usual constraints: 



G^v&v + G^X'tX* = , 

p,X'» = . 



(D.10) 
(D.ll) 



We then choose the usual gauge: p a = 2J = constant and t = r, where the total 
angular momentum in a is L = \J\ymJ ■ From (|D . 1 1| ) we can eliminate D\a and 
then use ( p,10| ) to eliminate pt from the Lagrangian (p.9| ) . The result is 



C = p^+p a X a + 2JC - VA ab PaPb + A 



(D.12) 
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where the indices a, b run trough the rest of the coordinates excluding t and a. The 
functions A ab and A are 



AJ 2 , r /2 _ 

/2 



A = t^ + 7^ + ^' 2 + (1 - r2 K , ( D -!3) 



2 



A " = r^ + i 4^ ! < 2 ' ( D - 16 ) 

A r0 = K 4>r = (1 ~ r2 YV ^ (D 17) 



(1 — r^r'n't 
Tj 2 
1 — r 2 )(j)'n', 

1j 2 
(1 - r^n^n'. 



A" = A ir = ^ ^ ' , (D.18) 

A^ = A* = {L ~ A '2 q} "' i > ( D - 19 ) 



A« = A" = 4J2 l J . (D.20) 

(D.21) 

We then proceed as before by eliminating the momenta using their equations of 
motion. This gives 

£ = 2JC - v/A(l - A ab x a x b ) , (D.22) 

where x a = X a for a ^ (p and x* = 1 + (ft, and A a & is the inverse of the symmetric 
matrix A ab . 

One can now expand as usual in the limit of large J assuming X a ~ 0(1/ ' J 2 \ 
To lowest non-trivial order we get the action (|6.3|) except that for the r dependence 
of the last term we get J d 2 a(l — r 2 ) 2 (dirii) 2 . 



E. Coherent states for the g-deformed algebra 

The g-deformed algebra has the following representation |62| . 



a \n) = \/[n +1] \n + 1) , a\n) = \/[n] \n — 1} , with [n] = . (E.l) 

Coherent states can be defined as eigenstates of the annihilation operator, and they 
are given by |fJ3| 

oo 



i) = {exp q (\z\ 2 ))- 1/2 J2^F^\n), with \z\<- , (E.2) 
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where [n]\ = [l].[2] . . . [n]. The g-exponential is defined by, 



exp q (x) 



^Ml = n- (i-gHi-^)' (E ' 3) 



As usual, these coherent states constitute a non-orthogonal and overcomplete basis. 
The overlapping between coherent states is 

/ I '\ exp (zz') 

(z\z ) = y/2 ■ ^ E - 4 ) 

( eX P 9 (k| 2 ) ex P 9 (kT)) 

The resolution of the identity takes the form, 



1 /" „ 1 f 2w /-V(l-9) 

2^ 



2 



d 2 q z\z) (z\ = — J dxj> J d q (r 2 )\z)(z\ = 1 . (E.5) 



where _D 9 is the open disk in the complex plane with radius 1/(1 — q), z — re 1 ^ and 
the integral over r is the so-called Jackson Integral, 

/ f{x)d q x = a{l-q)Y^q k f{q k a), (E.6) 

JO „_ n 



ra=0 

which satisfies J ( j Cl (exp (3 (x))~ 1 a;' 1 (i g a: = [n]\. 
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